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O . Abstract We study chromospheric oscihations including umbral flashes and 

1^ I running penumbral waves in a sunspot of active region (AR) using scanning 

1^ . spectroscopy in Ha and Ca ll 8542 A with the Fast Imaging Solar Spectrograph 

\ (FISS) at the 1.6 meter New Solar Telescope at Big Bear Solar Observatory. A 

bisector method is applied to spectral observations to construct chromospheric 
f^// Doppler velocity maps. Temporal sequence analysis of these shows enhanced 

high-frequency oscillations inside the sunspot umbra in both lines. Their peak 
frequency gradually decreases outward from the umbra. The oscillation power 
is found to be associated with magnetic-field strength and inclination, with 
different relationships in different frequency bands. 

o : 

^ , Keywords: Sunspots, Magnetic Fields; Oscillations, Solar; Chromosphere, Ac- 

c/^ ■ tive. 

1. Introduction 

> 

I Umbral flashes (UFs) and running penumbral waves (RPWs) are long-known 

sunspot oscillation phenomena that have been studied extensively. The first is 
more three- minute, the latter more five- minute in character. UFs were first 
discovered by Beckers and Tallant (1969) in Call H and K filtergrams and 
spectrograms from sunspots. It has been suggested that UFs brightness forms in 
. the local umbral gas during the compressional phase of a magneto-acoustic wave 

(e.^., Havnes, 1970; Rouppe van der Voort et al.^ 2003) and recently it has been 
confirmed by Bard and Carlsson (2010) from the NLTE simulations of Call H 
line. RPWs were first studied by Zirin and Stein (1972). They appear as a sort 
/\ . of continuation of some of the flash waves. Main difference between the UFs and 

the RPWs is that the UFs are best seen in Call H and K lines while RPWs 
are best seen in Ha Dopplergrams. Furthermore, UFs and RPWs have been 
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interpreted as different manifestations of the same oscillatory phenomenon in 
a sunspot, combining upward-shock propagation with coherent wave spreading 
over the entire spot (Zhugzhda and Dzhalilov, 1984; Rouppe van der Voort et a/., 
2003; Bloomfield, Lagg, and Solanki, 2007). The flashes are present only when the 
oscillation amplitude is sufficiently large (> 5 kms~^) but oscillatory motions 
are always present in nearly every umbra (Moore, 1981). 

However, it is still unclear what determines the characteristic properties of 
the oscillations in sunspot chromospheres. How do the waves, that are associated 
with these oscillations, propagate in the magnetized atmosphere? What is the 
spatial variation of frequency across the sunspot and in features of different 
spatial scales? The oscillations associated with waves are important to study 
since they carry some information on the properties of the region from which 
they originate and through which they propagate. Furthermore, the rapidly 
evolving field of local helioseismology, to investigate the sub-surface structure 
and dynamics, uses the oscillation properties. 

UFs and RPWs appear more clearly in Dopplergrams than in intensities. 
In many earlier studies the chromospheric oscillations were analyzed by con- 
structing Dopplergrams from filtergrams at fixed wavelengths in the red and 
blue wings. In this paper, simultaneous spectral observations in the Ha and 
Call 8542 A lines have been used to construct Doppler velocity maps (hereafter 
Doppler maps) from the bisectors of these two lines. The rest of the article 
is organized as follows: We describe the observational data in Section 2 and 
present the methods of analysis in Section 3. Results and discussions are given 
in Section 4. Finally, Section 5 is devoted to the summary and conclusions. 

2. Observations 

We observed AR NO A A 11242 on 30 June 2011. It was located at heliographic 
latitude 17° N and longitude 29° W, and consisted of a well-developed sunspot 
of southern polarity with mean magnetic field strength 1200 G. The sunspot 
was surrounded by a number of small magnetic fragments of northern polarity. 
AR was well observed by the FISS instrument (Chae et al, 2012), and the 
Helioseismic and Magnetic Imager (HMI: Schou et a/., 2012) on board the Solar 
Dynamics Observatory (SDO). 

FISS is a slit spectrograph taking spectra with rapid-scan capability. It ob- 
serves the solar chromosphere simultaneously in two spectral bands centered 
around the lines Ha and Call 8542 A (hereafter Call) with spectral resolutions 
of 0.045A and 0.064A, respectively. The pixel resolution at both wavelengths is 
^O.W pixel-^ 

The FISS data cube, with a field of view of 24^' of scan width (total number 
of scans 150) and 41^' of slit length, covered the entire sunspot of AR NOAA 
11242. The observations were taken relatively under good seeing conditions dur- 
ing 17:51:38- 18:55:42 UT, except for a gap during 18:04:45 - 18:18:10 UT, with 
the scan step sampling, timing, and cadence of 27 s, 130 ms and 30 s, respectively, 
in both spectral bands. 

In order to study the association between chromospheric oscillations and mag- 
netic fields, we have used the HMI vector-magnetograms: field strength, azimuth 
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angle, and angle of inclination. The angle of inclination is measured from the line- 
of-sight and azimuth from -\-Y direction of the heliocentric coordinates (X,Y) 
with from apparent disk center through the solar poles. 

3. Data Reduction 

Out of 105 cubes for Ha and Ca ll each, seven were discarded being of low quality 
from failure of adoptive optics locking. The data are processed in two main steps: 
pre-processing, and post-processing. In the first step, we corrected the data for 
bias, dark, flat and slit patterns. In the second step the wavelength scales were 
calibrated. Finally the processed data are compressed using principal component 
analysis (PC A: Pearson, 1901) which is useful specially when spectral profiles 
are similar to one another (Rees et a/., 2000). PC A compression suppresses the 
noise to a substantial amount without much loss of information. More detail on 
the processing and PC A compression can be found in Chae et al. (2012). 

3.1. Alignment 

Image rotation at the Coude focus, telescope guiding errors, errors in slit po- 
sitioning, and seeing cause shifts between successive scans. We used the far- 
wing parts of the scans as references. The successive scans were first derotated 
and then aligned by cross-correlation to the reference scan, and the Ha and 
Call scans were also aligned by cross-correlation. The error due to slit - posi- 
tioning is found to be ±0.16'^ Typical values of the net shifts in scan x-and 
^/-directions were found to be 8 pixel and 2 pixels, respectively. The FISS data 
were aligned with the HMI images, by manual feature and pattern matching, 
into solar (X, Y) coordinates. 

Figure 1 shows a mosaic of the aligned FISS and HMI images of AR NO A A 
11242 on 30 June 2011. The FISS Ha (a) and Call (b) images are constructed 
from the spectral data at wavelength 4. OA from line center. The bottom panel 
shows the HMI white light (c) and magnetogram (d) images of the same region 
observed at nearly the same time as the FISS images (specified in each panel 
bottom). The line PQ drawn in panel (a) marks the reference locations along 
which the characteristics of oscillations have been analyzed (Section 4). The 
arrow in the bottom- left panel shows the direction of solar disk center. 

3.2. Timing Correction 

Both the image rotation and the image shifts between successive scans that were 
corrected by derotation and alignment imply that the timing per pixel has no 
longer a fixed cadence. We determined the actual observing times per pixel from 
the alignment vectors. 

Figure 2 shows the temporal variation in cadence, (^(At), for 3 locations (1,2 
and 3 marked with circles in Figure lb). There is a data gap between 18:04:45 
- 18:18:10 UT, also the images have large shift and rotation. In addition, (5(At) 
varies in the range ±5s, for the total observation period of 63 minutes. The 
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Figure 1. Image mosaic of the AR NO A A 11242 on 30 June 2012 in different wavelengths. 
The upper panels show the raster images at wavelengths 4.0A from Ha (a) and Call (b) line 
center. The line PQ (a) marks reference locations for the analysis. The locations marked by the 
symbols "x" and "+" represent the approximate boundaries of the umbra and the penumbra, 
respectively. Small circles labeled with integer numbers (b) mark the locations for which time 
information is shown in Figure 2 while big circles labeled with A, B and C represent three 
locations as in Figure 4. The lower panels show the HMI white-light (c) and magnetogram (d) 
images observed near the time of the FISS observation (upper panels). The arrow in panel (c) 
denotes the direction of the solar disk center. 



deviation S{At) for the point 1 (3) becomes more negative (positive) with time 
which suggests that the pixel corresponding to this position was observed earher 
(latter). Note that we have only used the data sets starting from the time 
18:18:10 UT, for the oscillation studies. Also the 7 bad- AO scans were taken 
out. 
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Figure 2. Timing shifts due to image rotation and scanner shifts at three locations marked 
in Figure 1(b). Vertical dotted line specifies the reference image. Arrows in the upper panel 
mark the locations of bad- AO scans. 



3.3. Bisector Measurement, Doppler Maps and Power Maps 

We compute the Doppler velocity using the bisectors (Gray, 1976; Dravins, 
Lindegren, and Nordlund, 1981) of spectral line profiles. The locus of bisectors 
represents the asymmetry of the spectral line around the line center. For the 
bisector for a given spatial location, one may sample the profile asymmetry at a 
given intensity level (Cavallini et a/., 1986; Keil and Yackovich, 1981), at a given 
separation from the nominal line center wavelength or the per-pixel measured 
minimum wavelength (Bhatnagar, Livingston, and Harvey, 1972; von Uexkiill, 
Kneer, and Mattig, 1983; Tziotziou, Tsiropoula, and Mein, 2002; Tziotziou et a/., 
2007), and lambdameter measurement at given chord length (Slaughter and Wil- 
son, 1972; Stebbins and Goode, 1987). We prefer to measure Doppler velocities 
with the lambdameter method. 
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Of course lambdameter method has problem, along with the first two meth- 
ods, in the case when the Doppler width is changing, e.g.^ spectral lines are wider 
in the sunspot umbra than quiet region. There is a possibility of width variation 
with time for a given location due to energetic activities. However, it is difficult 
to remove the width related variation in the Doppler velocity for every pixel in 
the time series. To minimize this effect, we have chosen an optimum value of the 
bisector chord, SX = 0.35A(0.15A) for Ha(Call). These chords were applied to 
all the spatial pixels to compute the Doppler maps (see Figure 4, bottom row). To 
average over remaining image distortions due to seeing, guiding errors, etc., we 
applied spatial 3x3 boxcar smoothing. For the power analysis of Doppler maps, 
we prefer not to interpolate the data to equal temporal sampling but instead used 
the Lomb-Scargle periodogram (LSP) technique (Lomb, 1976; Scargle, 1982). 




Figure 3. Time- average intensity profiles (thick-solid) and bisectors (thick-dashed) for the 
location B (see Figure 4). Thin-dashed curves, plotted with respect to x-axes drawn near 
lines' cores, represent bisectors on a magnified scale (50 times for Ha and 20 times for Call) 
sampled at the same location at different times. The horizontal arrow marks bisector chord 
SX = 0.35A(0.15A) for Ha(Caii). 



4. Results and Discussions 

Figure 3 shows time-averaged intensity profiles (thick-solid) and bisectors (thick- 
dashed) of Ha (a) and Call (b) lines for location B (see Figure 1(b) and Figure 4). 
Thin-dashed curves show the magnified bisectors (50 and 20 times respectively 
for the Ha and Call) at different times for the same location B, and can be read 
from the x-axes drawn near the line cores. 

The bisectors for the Ha line is "C" shaped while for the Call line they 
are inverse of "C", and represent blue and red asymmetry, respectively, in line 
profiles. Uitenbroek (2006) has also reported inverse "C" shaped bisectors and a 
strong red asymmetry in the Call line. The amount of curvature or equivalently 
the velocity covered by the bisector is a measure of asymmetric flows. 
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Figure 4. Simultaneous intensity (upper row) and Doppler velocity (lower row) maps of AR 
NOAA 11242 computed using bisector method for the Ha((5A=0.35A) and Ca ll((5A=0.15A). 
Arrows mark the umbral locations where Doppler velocity measurement in Call failed due to 
emission features in umbral region. The positions A, B and C labeled with circles represent 
the locations for which wavelength-time maps are shown in Figure 5. 



Figure 4 (bottom row) shows samples of chromospheric Doppler maps of the 
Ha and Call lines respectively constructed at bisector chords = 0.35 A (c) and 
SX = 0.15 A (d), where the upper panels are intensity maps obtained from the 
mean of the intensities of the two wings at the chord intersections. The Doppler 
velocities in these maps range ±6kms~^, where negative/positive velocity shows 
upward/downward plasma motions. Patterns associated with running penumbral 
waves, in Ha and Call, can be seen in the Doppler maps. The overall red-blue 
asymmetry of the Doppler maps is caused by off-disk-center viewing. The blue in 
the limbward part of the penumbra indicates the line-of-sight component of the 
reverse Evershed flow. It is less evident on the center side, although filaments 
there would have better alignment because this spot has no penumbra in the 
centerside quadrant. 
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Near the umbra in Ca ll, there are pixels with high velocities (mark by arrow) 
which are artifacts and occur due to failure of the method adopted for the 
Doppler velocity determination. Similar kinds of Doppler velocity and magnetic 
field artifacts have also been reported earlier in the measurements obtained from 
the Michelson Doppler Imager (Maurya and Ambastha, 2009; Maurya, 2010) and 
the HMI (Maurya, Vemareddy, and Ambastha, 2012). We found that these pixels 
are associated with emission features of UFs marked by arrows in the intensity 
images. This problem was already described by Kneer, Mattig, and v. Uexkuell 
(1981), and also by Tziotziou et al (2006). UF emission features in Call can 
also be seen in Figure 5. 




-20 20 -20 20 -20 20 -20 20 -20 20 -20 20 
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Figure 5. Intensity as a function of wavelength and time for the Ho; and Call lines at three 
spatial locations, A, B and C, respectively for the umbra, penumbra and super-penumbra 
(see Figure 4). Dotted curves, passing through the line's centers, show the Doppler velocities 
(magnified by factor of 5 for all the panels except the panel 2 from the left) at corresponding 
locations. 

Figure 5 shows the intensity variations in Ha (panels, 1, 3 and 5 from the 
left) and Call (panels, 2, 4 and 6 from the left) as a function of time and wave- 
length for the three spatial locations A, B and C corresponding to the umbral, 
penumbral and outer-penumbral regions, respectively. The Doppler velocities 
for bisector chords lengths of 0.35A(Ha) and 0.15A(Call) for corresponding 
positions are over plotted (dotted curves). Here note that the computed Doppler 
velocities are much smaller than the spectral ranges as shown in the the abscissa 
axes, and therefore we magnified their values by factors of 5 for both the Ha and 
Call, respectively, except the panel 2 from the left. 

For the position A, the Ca ll line core shows emission features of UFs and the 
Doppler velocities computed through bisectors may be inaccurate. Nevertheless, 
the Doppler velocities for all locations vary quasi-sinusoidally with time which 
represents the oscillatory motions of the chromospheric plasma emitting the 
Ha and Call lines. There is an interesting pattern of different periods when we 
move from the panel A to C, corresponding to different locations. The period of 
oscillations increases from A to C. 
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Figure 6. Space-time maps of intensity (upper row) and Doppler velocity (lower row), in 
¥La (left column) and Call (right column), along the reference line PQ (see Figure la). Dotted 
horizontal lines around at 16^^ and 9^^ represent the approximate boundaries of the umbra and 
the penumbra, respectively. 



Figure 6 shows space-time intensity maps (upper row) for Ha-0.35A (a) and 
Call-0.15A (b), and Doppler maps (lower row) for the Ha(^A=0.35A) (c) and 
Caii((5A=0.15A) (d) along the line PQ (see Fi gure la). It is evident that the con- 
trast of the intensity and Doppler velocity variations in the umbral-penumbral 
region is large. Also the contrast in Ha is larger than in Call. In the umbral 
region, the ridges are almost vertical and show sinusoidal patterns of intensity 
and Doppler variations with time. Near the umbral boundary, there is a sharp 
change in the slope of the ridges. The intensity and Doppler velocity contrasts 
in the ridges decay with distance from the umbral boundary to outer edge of 
the penumbra, and it is hard to distinguish intensity ridges further. This shows 
that the RPWs decay in the super-penumbral regions of the sunspot. 

From the raster images, we found that the amplitude of RPWs gradually 
decreases from the line centers. It is difficult to identify the RPWs beyond the 
0.75A in the Ha and 0.50A in the Call. Also, the amplitude of RPWs decreases 
from the umbral boundary to outer edge of the penumbra which is consistent 
with the earlier reports (Giovanelli, 1972; Zirin and Stein, 1972). 
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Figure 7. Average oscillation power maps in different frequency bands in mHz (from the left 
to the right) for the iia{6X=0.35A) (upper) and Ca ll(5A=0.15A) (lower). 



The space-time Doppler panels of Figure 6 show that the RPWs have much 
more blue than red in Ha, the reverse in Call but no such asymmetry appears 
in Tziotziou et al. (2006, 2007). However, their traces also have larger amplitude 
for Call in the umbra plus a shift between largely negative for Ha and largely 
positive for Call. 

Figure 7 shows a sample of oscillation power maps, for the Ha(^A=0.35A) 
(upper) and Call((5A=0.15A)(lower), averaged over different frequency bands. 
For comparison, all the maps are shown in the same amplitude range (see color 
scales in the top). It is evident that there is a large power in the sunspot umbra 
at higher frequency band (4.5 < < 7.5 mHz) than the lower frequency band 
(1.5 < ly < 4.5 mHz) while in the outer penumbra the power is large at lower 
frequencies than in higher frequencies. This is consistent with earlier reports, 
e.g.^ Rouppe van der Voort et al. (2003). This pattern exists in both the spectral 
bands. The high frequency oscillations in Call are found to be concentrated in 
smaller regions than in Ha, although there are some bad data points in the umbra 
of Call maps due to artifacts in the Doppler velocity images {of. Figure 4). 

Figure 8 shows the oscillation power as a function of space (along the line 
PQ in Figure 1(a)) and frequency. There is a clear linear trend in the maximum 
power variation from the umbra to outward. The oscillations in the three and 
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Figure 8. Oscillation power maps along the line PQ (see Figure la) for the Doppler velocities 
of Ha((5A=0.35A) and Ca ll(5A=0.15A). Dotted lines mark the boundaries as shown in Figure 6. 



five minutes bands along with the features at other frequencies can be seen in 
both maps. 

In order to study the nature of the waves in the sunspot, we performed the 
Fourier phase difference analysis between the two Doppler signals obtained for 
Ha and Call. The Dopplershift measurements along PQ were interpolated to 
fixed- interval (30 s) sampling for this purpose. 

Figure 9 shows the phase difference between the Ha((5A=0.35A) and Call 
((5A=0.15A) Doppler velocities as function of frequency for four locations in the 
umbra(a), penumbra(b) and superpenumbra (c); locations are marked in the 
panel (a) of Figure 10 by symbols, □, O, x, and +. Symbol sizes represent the 
cross-spectral power between the two Doppler velocities. The average phase dif- 
ference in the frequency bands 5 and 3 mHz are around 40° where cross-spectral 
powers are significant. The phase difference shows that the UFs and RPWs 
are caused due to upward propagating MHD waves. The phase difference in the 
superpenumbral region is about 10° at lower frequency ~1 mHz where significant 
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Figure 9. Phase difference (Acp) in Ca ii((5A=0.15A) and iia{SX=0.35A) Doppler velocities 
of: (a) umbra, (b) penumbra, and (c) superpenumbra, at four locations marked by symbols, 
□ , O, X, and + along the reference line PQ (Figure 10a), where Call((5A=0.15A) is taken as 
leading. The distance d is measured from the point P along PQ. The symbol sizes represent 
the cross-spectral power between the two Doppler time series. 



power exists. This difference may caused due noise and unequal sampling in the 
original data. 

In order to ascertain the spatial variation in the oscillation power from umbra 
to outward, we have computed the root-mean-square (RMS) velocity (vrms) along 
the line PQ and results are shown in Figure 10. Dotted vertical lines correspond 
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Figure 10. Variation m frms 

velocities of Ha((5A=0.35A) and Caii((5A=0.15A) along the 
reference lines QN(b), QM(d) and QL(c) marked in the panel (a). Vertical dotted lines show the 
approximate umbral and penumbral boundaries marked by symbols "x" and "+", respectively, 
over the solid curves in the panel (a). Other symbols "□C>x+" in the umbra, penumbra and 
superpenumbra of the sunspot mark the locations along the line PQ (same as in Figure la) 
for phase analysis in Figure 9. 



to approximate umbral and penumbral boundaries marked by symbols "x" and 
respectively, in the panel (a) of Figure 10. It is evident that Vrms is large 
in the umbra of both the Ha((5A=0.35A) and Call((5A=0.15A) and decreases 
(increases) in the penumbra (superpenumbra). In the inner umbra, Call Vrms for 
point Q rises considerably above the 2.1kms~^ axis cutoff which is caused by 
combined effects of high Doppler velocities, and large uncertainty in the Doppler- 
shift measurements. In the outer umbra, Vrms of Ha is larger than Vrms of Call. 
This may be due to the difference in formation height between the two spectral 
lines or other difference in response. For example. Call may have more opacity 
in cool post-shock gas than Ha (which gets large opacity only in hot gas), so 
that Call feels cool post-shock downdrafts and Ha does not. Call then may get 
the K2V-like shock grain pattern from shock interference along the line of sight. 
Towards P the LOS alignment is lost. The large Vrms in umbra consists mostly of 
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high-frequency power. The high-frequency components hkely describe sawtooth 
shocks (see Figure 5). 



o 



1.5^i/<2.5 





1.5^u<2.5 •• ■ 











?^ 

CD 
O 



2.5^i/<4.5 





2.5gi/<4.5 : 















6 

^ 4 

CD 




4.5^i/<6.5 



— 6 

^ 4 
?^ 

CD 



• Hc(((5X = 0.35A) 

A CaII((5A = 0.15A) 






6.5^i/<7.5 - 






- ' • 








20 40 


60 80 


7 (de^ 


^ree) 



5 10 15 

B (xlO^ G) 

Figure 11. Variation of the oscillation power in different frequency bands (from top to the 
bottom) with: (left) total magnetic field (B), and (right) inclination angle (7) measured from 
a rectangular area along line PQ of width 10 pixels. 
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Figure 11 (left column) shows the association between oscillation power and 
corresponding magnetic field strength B from a rectangular area along line PQ 
of width 10 pixels. At lower frequencies (1.5 < v < 2.5 mHz) the power is largest 
in areas with smaller B which correspond to the regions of outward penumbra, 
and rapidly decreases with B. In the frequency band 2.5 < u < 4.5 mHz, the 
power initially increases up to B ^ 700 G and then decreases with B. In the 
frequency band 4.5 < f < 6.5 mHz the oscillation power increases with B and 
becomes maximum around 1600 G then decreases rapidly. In the frequency band 
6.5 < V < 7.5 mHz most of the power is concentrated near the very high B 
(> 1500 G) which correspond to the umbral region of the sunspot while the 
power is nearly constant for the B range 300-1200 G. 

For the study of association between the oscillations and the field inclinations, 
we computed the inclination from the vertical to the surface using relation 
7 = cos~^(5r/|B|), where, Br is the radial component of the magnetic field B. 
Uniform shear method is used (Moon et al.^ 2003) to resolve the 180° ambiguity 
in the azimuth angle. 

Figure 11 shows the variation in the oscillation power with magnetic field 
inclination 7 corresponding to the left column of magnetic field strength B. It 
is evident that the power shows reverse trend with inclination angle than the 
field strength. In the lower frequency band (1.5 < v < 2.5 mHz), most power is 
concentrated at higher 7 > 80°. The peak of power distribution shifts toward 
lower inclination for higher frequency bands. For instance in the frequency band 
6.5 < ly < 7.5 mHz, most of the power is concentrated in the inclination region 
10°- 30°. 

The opposite relation of inclination and field strength with the power is 
obvious from the magnetic field distribution in the sunspot umbra; we have 
strong field with small inclination. From the umbra to outward the field strength 
increases while inclination decreases. Figure 11 shows that the oscillation power 
is associated with the field strength and the inclination angle both but they 
behave differently in different frequency bands. 

5. Summary and Conclusions 

We studied the nature of running penumbral waves using Ha and Ca ll Doppler 
images constructed with a bisector method from area scans with the FISS 
spectrometer. We found that the RPWs are easily seen in the intensity images 
constructed near the core and their amplitude decay with wavelength from the 
line center to outward. Also their amplitude decreases from the umbra to out- 
ward. These results are consistent with earlier results about the RPWs (Zirin 
and Stein, 1972; Christopoulou, Georgakilas, and Koutchmy, 2000; Tziotziou 
et a/., 2006). The running penumbral waves decay near the outer boundary of 
penumbra. 

We found that the chromosphere umbra of the sunspot shows large RMS 
velocity ('^rms) in both spectral lines. From the umbra to outward the Vj-ms grad- 
ually decreases, and there is no distinct boundary between the v^ms of the umbra 
and penumbra, however, Vj-ms has minimum in the penumbral region. There is 
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another interesting pattern seen in the Vj-ms of Ha and Call hnes, Vj-ms is smaher 
in the Call than Ha away from the umbral center. These results reveal that the 
Ha and Call are formed differently in the shocks that make up UFs. Furthermore 
Vrms is related to the oscillation power (power ex '^rms^)- That is, we have larger 
power in the umbra than in penumbra, also, the power is larger in the Ha than 
in Call. 

Our times series analysis of chromospheric Doppler maps shows high - fre- 
quency power in the umbra of both the spectral lines. Ha and Call, which is 
also evident from the Vrms distribution. This study shows that the total power 
gradually decreases with frequency from the umbra to outward which confirms 
earlier reports (Tziotziou et al, 2006; Nagashima et aL, 2007; Socas-Navarro 
et al.^ 2009). The power maps indicate that this decrease is set by the high- 
frequency contribution. Our analysis of spectral observations shows that the high 
frequency oscillations exist in the umbra of 5 min band. The observed power in 
the penumbral region shows strong peak in the 3 min band while the power is 
small at other frequencies. 

The penumbral oscillations are regarded as a tail of five minutes oscillations 
resulting further dependence of acoustic cutoff frequency on the magnetic field 
inclination (Cally, Bogdan, and Zweibel, 1994; De Pontieu, Erdelyi, and James, 
2004). The oscillation power shows a strong relationship with the magnetic field 
strength and angle of inclination. Our analysis showed that the peak oscillation 
frequency not only depend upon the inclination but also on the field strength, 
although they behave differently in different frequency bands. 

Observations in Ha and Call bands demonstrated the oscillations with fre- 
quencies around 4.8 mHz are dominant in the umbral and inner penumbral 
regions, together with the presence of RPWs in the band of 3 mHz. Although 
the scale height between chromospheric Ha and Call lines are closer to the 
photosphere, the relationship of frequency and magnetic fields may be affected. 

We can summarize the above results as follow: 

• Call showed reverse "C" shaped bisectors at some locations of the sunspot. 

• Bisector- measured Doppler velocity in Ca ll is affected by umbral emission 
features. 

• Amplitude of RPWs decreases with distance from the umbra to outward 
and decay near the penumbral boundary. 

• RPWs have much more blue than red in Ha and reverse in Call. 

• Peak power frequency gradually increases from umbra to outward in both 
lines. 

• In umbra and penumbra of sunspot. Ha and Call lines showed a phase 
difference of ~ 40° in the 3- and 5-minute oscillation bands. 

• Vrms in the umbra is larger than penumbra of both lines. In outer umbra, 
Vrms is larger in Ha than Vj-ms in Call. 

• The oscillation power shows different relation with the magnetic field strength^ 
and the inclination in different frequency bands. 

We conclude that the main oscillations properties still suggest MHD waves 
propagating upward along fanning field, as in earlier studies, but with inter- 
esting but yet unexplained response differences between the two lines. The 



SOLA: Chrom0sc_v2.tex; 7 February 2013; 2:25; p. 16 



Chromospheric Sunspot Oscillations in Ho; and Call 8542A 



transformation of vertical oscillations into RPWs, at the boundary of umbra 
and penumbra, is an interesting feature which require further observations and 
theoretical studies to explain the wave phenomenon in the sunspot. 
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